Excitation Spectra of the Negative-f/ Hubbard Model: 
A Small- Cluster Study 
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An exact-diagonalization technique on small clusters is used to study low-lying excitations 
and superconductivity in the two-dimensional negative- U Hubbard model. We first calculate the 
Bogoliubov-quasiparticle spectrum, condensation amplitude, and coherence length as functions of 
the coupling strength U/t, thereby working out how the picture of Bogoliubov quasiparticles in 
the BCS superconductors is affected by increasing the attraction. We then define the Cooper-pair 
operator as a spatially-extended composite boson and make a variational evaluation of its internal 
structure. We thereby calculate the single-particle spectral function of the Cooper pair and obtain 
the dispersion relation for its translational motion. The dynamical density correlation function of 
the pairs is also calculated. We thus demonstrate the applicability of our numerical method for 
gaining insight into low-lying excitations of models for the intermediate coupling superconductivity 
relevant to cuprate materials. 
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I. INTRODUCTION 

High-temperature superconductors typically show very 
short coherence length and small carrier numbers, and 
may be in an intermediate regime between two well- 
understood regimes, i.e., the BCS weak-coupling su- 
perconductivity and Bose-condensation of pre-formed 
bosons. Although much theoretical effort has been de- 
voted to examining this crossover regime, not very much 
is known, in particular, of the low-lying excitations [1,2]. 
To consider this problem, we have studied the negative- 
U Hubbard model because of its controllable strength of 
attraction. This model, much discussed recently [3-10], 
may be a good reference system to more relevant models, 
such as the t—J model, for cuprate superconductivity. 

The negative-?/ Hubbard Hamiltonian is written 



H = —t 



E (4 



H.C. 



n^n n 



(1) 



<lj>cr 



with a negative value of U, where C; (cv) is the elec- 
tron creation (annihilation) operator at site i and spin a 
(=t,|), and ni a =c\ a Ci a . The summation <ij> is taken 
over all the rearest-neighbor pairs on the two-dimensional 
square lattice. We employ an exact-diagonalization tech- 
nique on small clusters; a 4x4 cluster with periodic 
boundary condition is used throughout the work, and 
ground states and dynamical correlation functions are 
calculated by the Lanczos algorithm. 

As an approach from weak-attraction limit we first cal- 
culate the Bogoliubov quasiparticle spectra in the BCS 



pairing theory as a function of U/t, and demonstrate how 
the picture of Bogoliubov quasiparticles is affected by in- 
crease of the strength of attraction [11]. Condensation 
amplitude and coherence length are also estimated. We 
then make a variational calculation of the Cooper-pair 
wave function and examine U ft dependence of the inter- 
nal structure of the pair. Next we calculate the wave 
functions of pairs with non-zero momenta and examine 
the single-particle excitation spectra of the pair to see 
its translational motion. Finally we calculate the dy- 
namical pair-density correlation, and examine a density 
fluctuation of the pairs. We thus demonstrate the ap- 
plicability of our numerical method for gaining insight 
into low-lying excitations of models for the intermediate 
coupling superconductivity. 



II. BOGOLIUBOV QUASIPARTICLE 

Let us first examine the validity of the Bogoliubov- 
quasiparticle picture in the negative-?/ Hubbard model as 
a function of attraction U/t. We use a recently proposed 
technique [12], i.e., an exact calculation of Bogoliubov 
quasiparticle spectrum on small clusters. We thereby see 
whether low-lying states of the negative-f/ Hubbard clus- 
ter of a given attraction U/t can be described by this 
picture. 

We define the one-particle anomalous Green's function 

as 
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where \ipg} is the cluster ground state with N (even) 
electrons, E is the ground-state energy (averaged over 
energies of the N and iV+2 electron states), and cj^ is 
the Fourier transform of c\ . We then define the spectral 
function F(k, oj)=— (l/w)lmG(k, oj+irj) with i]=+0 and 
its frequency integral 
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This Green's function should describe the excitation 
of a Bogoliubov quasiparticlc in the cluster, i.e., 
F(k, ijj)~F-^5(uj— Fk) with F k =Ak/2Fk for the quasipar- 
ticle energy Fk and gap function Ak, provided that the 
low-lying states of the cluster can be described by the 
microcanonical version of the BCS pairing theory [13]. 

The calculated results for F(k, u) at N~8 and 10 are 
shown in Figs. 1 (a) and (b), respectively. Fermi mo- 
mentum k F is at (7r/2,0) for N=8 and at (ir/2,n/2) 
and (?r,0) for N=10. We find that, at U/t=-2, the 
spectrum fits fairly well with the BCS spectrum (dotted 
curves) with s-wave gap function Ak=A : pronounced 
low-energy peaks appear at kF, smaller peaks appear at 
higher energies for other momenta. The observed gap 
is isotropic, and the spectral weights are consistent with 
the BCS form of the condensation amplitude Fk with a 
maximum at kp (see Fig. 2). However, at U/t=—6, the 
peaks already extend to higher-energy regions and the 
dispersion becomes rather deformed, indicating that the 
Bogoliubov quasiparticlc is becoming a less well-defined 
excitation although the lowest-energy peaks at kp are 
still sharp and well-defined. For much larger values of 
\U\/t the spectra are totally incoherent and the notion of 
Bogoliubov quasiparticles loses its significance. 

We also calculate the coherence length £ defined as 



e 2 = Ei v ^i 2 /Ei^ 
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The result (see Fig. 2) shows a rapid but smooth 
crossover from a Cooper-paired state to a Bose condensed 
state of tightly bound pairs. At \U\/tc±3— 4, £ is of the 
size of the nearest-neighbor to next-nearest-neighbor dis- 
tance. Note that Fk goes to zero when \U\/t— >0 for k^kp 
in agreement with the BCS pairing theory; nondivergence 
of £ at \U\/t—*Q is an obvious finite-size effect. 



III. COOPER PAIR 

As an approach from strong-attraction regime we in- 
troduce a spatially-extended 'composite-boson' operator 
defined as b\ = J2 r a ( r ) c iT c i+ri [14] and its Fourier trans- 
form & q with a(k) = ^ r a(r)e lk r . We then define a vari- 
ational wave function by adding a zero-momentum pair 
of up and down electrons to the iV-electron ground state: 
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where N s is the number of lattice sites. We determine 
variational coefficients a(k) so as to maximize the over- 
lap with the exact wavefunction \{ipQ +2 \4>o +2 )\ 2 under 
the condition that +2 ) be normalized to unity, which 
leads to a generalized eigenvalue problem 



Y}.KFw - AAWMk') = 
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with F k of Eq. 3 and N^ k i={^\c ki C-\ s .^_ v ^, ] \ [ ^) 
[15]. We find for the cluster that an overlap of A~0.95 
is generally achieved. Defining Fk by replacing tpo^ 
with iJ)q +2 in Eq. 3, we have Fk=Ek' a * (k')-^kk' and 
Ek -Fka(k)=l. One may use these relations with the 
exact Fk to check the validity of the variational wave 
function; actually we find values of ^ k Fk«(k) (=A X / 2 ) 
to be close to 1 within a few percent. 

Figure 3 shows the Cooper-pair wave function thus de- 
termined; the pair is added to the state with 7V=8 and 
k F =(7r/2,0). We find in Fig. 3 (a) that when |U|/i<6 
a(k) has a peak at k F , and with increasing \U\/t, it 
broadens over the entire Brillouin zone; the Pauli prin- 
ciple acting between the added electrons and Fermi sea 
plays an essential role in the structure of the Cooper-pair 
wave function. In real space, as is expected from the BCS 
pairing theory, a(r) (see Fig. 3 (b)) shows an oscillation 
of wavelength 2tt /fop corresponding to the peak in a(k), 
and the oscillation decays with the length scale of £ as 
one may see via U/t dependence of a(r). 

IV. EXCITATIONS OF COMPOSITE BOSONS 

Let us next examine excitation spectra for transla- 
tional motion of the Copper pair. We first introduce a 
creation operator of the finite-momentum (q) pair which 
is defined as 
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with q-dependent internal structure a q (k) (which re- 
duces to 6 q defined in Sec. Ill at q=0), and again perform 
the variational evaluation of a q (k) by a generalization of 
the above procedure where Fk and 7V k k' also depend on q. 
The operator thus obtained is compared below with the 
on-site pair-field operator 6 J. = /TvT^i V. A. A p- 1 ^ 1 



i.e., Eq. (7) with a q (k)=l. We find for the cluster with 
7V=8 that the overlap A>0.7 (mostly around ~0.8) is 
achieved by the optimization; the presence of the Fermi 
sea is again found to play an important role in the k and 
q dependence of a q (k). 

We then define the pair-addition spectrum as 
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whereby examining the single-particle excitation of the 
Cooper pair. The chemical potential for the Cooper 
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pair is defined as ^i p —dE/dN p —E +2 — Eq where N p 
(=N/2) is the number of pairs. Figure 4 shows the cal- 
culated result for P(q, ui) with optimized values of a q (k); 
we compare this with the result obtained for the on-site 
pair-field operator. We first of all find that with the op- 
timization higher-energy peaks are suppressed strongly 
and the lowest-energy peak is enhanced at each momen- 
tum q, indicating that the spatially-extended 'composite 
boson' describes the low-energy excitation of the system 
fairly well [16]; its effective mass, e.g., is noted to become 
heavier with increasing \U\/t. An indication of the off- 
diagonal long-range order is also seen in the calculated 
momentum distribution (6 q 6 q ), i.e., its enhancement at 
q=(0, 0), although in the thermodynamic limit a macro- 
scopic number of bosons are condensed into q=0 state. 
One may then check how the Bogoliubov theory for inter- 
acting bosons [17] works for these composite operators; 
we directly calculate the density correlation function de- 
fined as 



were partly carried out in the Computer Center of In- 
stitute for Molecular Science, Okazaki National Research 
Institutes. 
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with the density operator d q =(l/y / ]V7) ^ k 6^ 6k- The 
calculated result for B(q,ui) at U/t=—4 is shown in 
Fig. 5; the spectrum is approximated well by a single 
peak around q=(0,0), which is consistent with the ex- 
pected picture of the collective sound mode. For larger 
momenta this mode decays into pair-breaking excitations 
observed in F(k, ui). 



V. SUMMARY 

We have proposed numerical techniques for exam- 
ining the low-lying excitation spectra of the two- 
dimensional negative-?/ Hubbard model and studied the 
(i) Bogoliubov-quasiparticle excitation, (ii) condensation 
amplitude and coherence length, (iii) internal structure 
of the Cooper pairs, and (iv) single-particle and density 
excitations of the pairs, as functions of the attractive 
interaction. We have demonstrated how the picture of 
Bogoliubov-quasiparticle excitations in the BCS super- 
conductors loses its significance as the coupling strength 
increases and how the picture of spatially-extended 
composite bosons works as the coupling strength de- 
creases. Application of the present numerical technique 
to strongly-correlated electron models will help us gain 
insight into the intermediate-coupling superconductivity 
relevant to cuprate materials. 
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FIG. 1. Bogoliubov quasiparticle spectra F(k,a>) for the 
4x4 cluster with filling of (a) N=8 at U/t=-2 (left panel) 
and -6 (right panel), and (b) iV=10 at U/t=-2 (left panel) 
and —6 (right panel). Dotted curves in the left panels show 
the BCS spectral function obtained for (a) Ao/i=0.28 and (b) 
0.4. We use the value r)/t=Q.15 for Lorentzian broadening of 
the spectra. 

FIG. 2. Condensation amplitude F k and coherence length 
£ as a function of the attraction U/t calculated for the 4x4 
cluster with filling of N—8. 

FIG. 3. (a) Cooper-pair wave function a(k) and (b) its 
Fourier transform a(r) calculated for the 4x4 cluster with 
filling of N=8. 

FIG. 4. Pair-addition spectra P(q, u>) calculated for the 
4x4 cluster with filling of iV=8 at U/t=—4. Vertical dot- 
ted lines indicate fi p . Left panel: the spectra for b\ with 
optimized a q (k) values. Right panel: the spectra for the 
on-site pair-field operator. Lorentzian broadening of rj/t=0.05 
is used. 

FIG. 5. Density correlation function B(q, u>) of composite 
bosons calculated for the 4x4 cluster with iV=4 and U /t=—4. 
Lorentzian broadening of j)/t=0.05 is used. 
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